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ABSTRACT 

We investigate the success and problems of MOdified Newtonian Dynamics (MOND) in explaining the 
structural parameters and dynamics of remote Galactic globular clusters (GCs) and dwarf spheroidal 
(dSph) galaxies. Using the MOND value for the mass of the Milky Way as inferred from the Galactic 
rotation curve, we derive the tidal radii of Galactic GCs, and compare to observed values. Except for 
Pal 14, the predicted tidal radii of GCs are systematically larger than the observed nominal values. 
However, the associated uncertainties are so large that tidal radii are consistent on the la level. We have 
considered the importance of the Galactic tidal force on the survival of dSphs under MOND. Assuming 
mass-to- light ratios compatible with a naked stellar population, we found that the present Galactic dSphs 
preserve their integrity over one Hubble time, except Sextans which may survive the tidal interaction 
only for several Gyr. Mass-to-light ratios as inferred from the internal kinematics of dSph galaxies can 
be used, but they appear too large to be accounted for only by the stellar population in Willman 1, 
Coma Berenice, Ursa Minor, Draco, Ursa Major and possibly Bootes dwarves. Finally, the ability of 
the Sculptor dwarf to retain the observed population of low-mass X-ray binaries (LMXBs) is examined. 
Under the MOND paradigm, we find that the retention fraction in Sculptor is likely not larger than a few 
percent. Compared to the dark matter scenario, it turns out that MOND makes the retention problem 
worse. We propose that measurements of the radial velocities of the observed LMXBs in Sculptor could 
provide a way to distinguish between modified gravities or extended and massive dark matter halos. 

Subject headings: galaxies: dwarf - galaxies: individual (Sculptor) - galaxies: kinematics and 
dynamics - gravitation 



1. INTRODUCTION 

The dynamical mass-to-light ratios, M/L, derived in 
galaxies are usually larger than the expected mass-to-light 
ratio of the stellar component, usually interpreted as indi- 
cating that they must contain dark matter. Alternatively, 
one could argue that the discrepancy between total mass 
and baryonic mass is telling us that the Newtonian law 
of gravity is not governing the dynamics. In particular, 
the Mod i fied N ewtonian Dynamics (MOND) proposed by 
iMilgroinl (jl983[ ) has been proven to be successful in ex- 
plaining how galaxies rotate, without any dark matter (see 
Sanders & McGaugh 2002 for a review) . It is important to 
note that MOND ideas and variants have so far been intro- 
duced, tailored and calibrated with the explicit objective 
of accounting for the first order gravitational effects of hy- 
pothetical dark mater halos: rotation curves in large galax- 
ies and stellar velocity dispersions in spheroidal systems. 
One should also note that a wide range of other gravita- 
tional effects appear in galactic dynamics, which in prin- 
ciple, offer independent restrictions on any proposed mod- 
ified dynamic scheme. Some examples include the grav- 
itational stability of galactic disks (e.g., Sanchez-Salcedo 
& Hidalgo- Gamez 1999) or relaxat ion processes in stel- 
lar systems (jCiotti fc BinnevI l2004f ). Concerning dSphs, 
in Sanchez-Salcedo et al. (2006) we investigated the ef- 
fects the enhanced gravitational relevance of the stellar 
population under MOND would have upon the problem 
of the orbital decay of globular clusters (GCs) due to dy- 
namical friction, and found that although such a proposal 
can account for the velocity dispersion measurements, in- 



s piraling times become unco mfor tably shor t. 

iGerhard fc Spergell (|1992[ ) and iGerhardI (|1994) derived 
the mass-to-light ratios in M OND for s even dwarfs from 
virial arguments. Later on, iMilgronJ ([1995f) reanalyzed 
the M/L values with new data. He concluded that when 
uncertainties in the observed luminosity and central ve- 
locity dispersion are included, the values agree with those 
expected for a normal old stellar population. The flux 
of data about the kinematics of Local dSphs is growing 
rapidly. The new determinations of line-of-sight velocity 
profiles, the recent discovery of a dozen dSph galaxies, 
and the detection of low-mass X-ray binaries (LMXBs) in 
Sculptor may give crucial clues on the dark matter prob- 
lem and its alternatives. In the Newtonian approach, these 
observations lend support to the cosmological-motivated 
interpretation that some dSphs reside within the most 
massive and extended substru cture dark ma tter halos 
(iKlevna et al.l f2001. 2005; Dchn en fc KinEjl2006D . The fact 
that some dSph galaxies possess two distinct populations, 
as inferred from studies o f resolved stellar populations, 
e.g. lHernandez et al.l (|2000t ) for Carina and Leo I, also sug- 
gests an extended dark matter halo in order to prevent the 
loss of all the gas during the first starburst. The question 
that arises is wh ether MOND can mimic an extended halo 
in dSphs or not. iLokasI (|2002[ ) and lLokas et~all (|2006D have 
noticed that MOND experiences a serious difficulty in ex- 
plaining the velocity dispersion profile of Draco without 
dark matter. 

A wide view of the successes and problems of MOND 
in dSph galaxies in comparison to the dark matter sce- 
nario can motivate new tests on the nature of dark matter 



1 



2 



and its alternatives. Here we will investigate the ability 
of MOND in explaining the observed tidal radii in Galac- 
tic GCs and dSph galaxies (§2), and the problem of the 
retention of LMXBs in dSphs (§3), within the alternative 
optics of MOND without dark matter. As a byproduct, an 
updated compilation of the mass-to-light ratios in satellite 
dSphs under MOND is given. 

2. TIDAL RADII OF GCS AND THE SURVIVAL OF DSPHS 

In any law of gravity, the global M / L ratio in a spherical 
system can be derived in different ways; either modeling 
the observed velocity dispersion profile of the stars or us- 
ing the observed King tidal radius. The observed tidal 
radii can provide lower limits for the M / L that are inde- 
pendent of kinematic data (e.g., Faber & Lin 1983; Oh, 
Lin & Aarseth 1995; Pryor 1996; Burkert 1997; Walcher 
et al. 2003). 

The dynamical iKingl (Il962l ) tidal radius of a GC or a 
dSph orbiting in a circular orbit of radius D around a 
galaxy with a logarithmic potential is defined as: 

where Mdw is the mass of the satellite, Mg{D) the host 
galaxy mass inside D and k is Keenan's (1981a, b) factor, 
which depends on the gravity law, and accounts for the 
elongation of the zero-velocity surface along the line join- 
ing the cluster and galactic centers. For eccentric orbits. 
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where e and a are the orbital eccentricity and the semi- 
major axis, respectively. This formula for the tidal ra- 
dius is valid for both Newtonia n (e.g., Oh, Lin fc Aarseth 
1992) and MOND gravities (iGerhard fc Spergell 11992 
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provided that the masses and the Keenan factor are taken 
appropriately. In the MONDian case without dark mat- 
ter, Mdw and Mq are the "true" baryonic masses of the 
dSph satellite and the host galaxy, respectively. Knowing 
the Galactic mass radial profile and the eccentricity e and 
by equating King's tidal radius, rt, with the observed lim- 
iting radius, a lower limit on M/L of the satellite system 
can be inferred. 

The nominal value A: = 1, as first suggested in King's 
paper (1962), corresponding to the distance of the last 
closed zero- velocity surface along the line joining the cen- 
ter of the galaxy with the center of the satellite system, is 
commonly adopted as a useful fiducia l length scale under 
Newto nian gravity, after the work o f lOh. Lin fc Aarse"thl 
(|1992[ ) and lOh. Lin fc AarsethI (|l995t) . who demonstrated 
that even in the absence of significant two-body relaxation 
processes, the density profile of GCs in circular orbit re- 
sembles the King-Michie model with an asymptotic lim- 
iting radius consistent with the tidal radius as defined by 
[King (1962), with k^l. If wc adopt the same definition of 
tidal radius in MOND, i.e. the distance to the zero- velocity 
surface along t he axis pointing in the radial direction, we 
also get A; = 1 (jZhao fc Tianll2006[ ). 

The parameter F, defined as the ratio between the ob- 
served cutoff radius and the dynamical tidal radius as 



defined in Eq. ^ with k = 1, is a useful indicator of 
the importance of the Galactic tidal force on the inter- 
nal structure of GCs a nd dSphs (jOh. Lin fc Aarsethlll995l : 
iPiatek fc Prvoij|1995l ). For values F < 1, the tidal force 
is unable to stretch the satellite object significantly, and 
only internal processes may be considered. In this case, the 
stellar component of the dSph galaxies is not tidally trun- 
cated, and mass estimators based on the tidal radius un- 
derestimate the total M/L substantially. At values F ^ 2, 
a dSph can survive the tidal interaction for several Gyr. 
At these values a dSph could actually be unbound but not 
yet dispersed, as the timescales for disruption exceed sev- 
eral Gyrs. With the effect of the force increasing rapidly 
with increasing F, for _F ~ 3 the satellite disintegrates in 
a few orbits. Since it is unlikely that we are observing 
satellites in such a short-lived evolutionary state (except 
Sagittarius), F < 2 is expected. 

The values of F under MOND, Fm, for Galactic GCs 
and sa tel lite dSphs ha ve been calculated recently bv lZhaol 
(|2005D . IZhaol (|2005D uses a Keenan factor k = V2/3 
and equates the mass ratio with the luminosity ratio, 
i.e. Mdw/Mc « Ldw/Lc in Eq. (P), where Ldw and Lq are 
the luminosities of the dwarf galaxy and our Galaxy, re- 
spectively. Th e large scatt er found in the values of Fm was 
interpreted bv lZhaol (j2005D as challenging for any baryonic 
MOND Universe. 

In this Section, we revisit the consistency of the predic- 
tions of MOND with observations using the tidal radius 
theory, assuming the mass of the Milky Way to be known. 
1/3 We discuss GCs ( W2.1\i and dSphs ( W2.2^ separately, since 
remote GCs are likely tidally truncated, whereas dSph 
galaxies are not necessarily tidally truncated, even though 
their luminosity profiles can be fitted by King models. 

2.1. Globular clusters 

In order to estimate the tidal radii of Galactic globu- 
lar clusters under MOND, we must start by adopting a 
mass mo del for the Milky Way un der MOND. We take 
this from iFamaev fc BinnevI (j2005D (their model I) who 
use the Galactic rotation curve, over the region where it 
is well measured, the first 10 kpc, to derive a total bary- 
onic mass of Mg — 0.6-0.8 x 10^^ Mq (see also Dehnen fc 
Binney 1998). This yields a circular rotation at the solar 
circle of 200 km s~^ and good agreement with dynami- 
cal measurements at all radii internal to 10 kpc. Since 
the asymptotic circular velocity of test particles around a 
mass Mq is V = {GMcaoY^^, this imphes an asymptotic 
circular velocity of 170 ± 5 km s~^ (essentially shared by 
most of their models, see also Famaey et al. 2007), for our 
adopted qq = 0.9 x 10~^ cm s~^, which is the value derived 
in Begeman et al. (19 91) rescaled to the new distance scale 
(jBottema et al.|[2002f) . This mass compares quite well with 
the expected baryonic disk mass for a galaxy with a disk 
V band luminosity of ^ 1.5 x 10^° L©,^- The asymptotic 
value of the rotation curve might seem a little on the low 
side, but one must bear in mind that more familiar values 
of order of 200 km are the result of extrapolating dark 
halo mass models beyond the regions where the rotation 
curve as such is actually measured. 

It is worthwhile to compare the King tidal radii in the 
Newtonian dark matter scenario and in MOND. According 
to Eq. ([T]), the tidal radius in the conventional dark mat- 
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ter model depends on the total mass of the host galaxy as 
[MGMr{D)+MG,dm{D)Y/'^, where Mgmt and Ma Am are 
the baryonic and dark mass, respectively. In MOND, the 
tidal radius goes as [Mbar{D)Y^^ ^ where M^ar is the bary- 
onic mass only. So that, remote GCs (those at D > 35 
kpc) should be the best discriminators between MOND 
and Newtonian dynamics as the dark halo mass becomes 
significant at large distances. 

Unfortunately, the range Z? > 35 kpc is mainly popu- 
lated by sparse, low-luminosity clusters with rather un- 
certain tidal radii and integrated magnitudes. For il- 
lustration, for NGC 2419, which is the fourth most lu- 
minous cluster of the Galaxy, discrepant values for the 
global Ml L and cutoff radius can be found in the liter- 
ature. iPrvor fc MevlanI (| 19931 ) report M/L = 1.2 and 
n = 230 pc (see also Trager et al. 1995), whereas Ol- 
szewski et al. (1993) derive a global M/L rati o of 0.7 ±0.4 
and T t = 280 pc. As discus sed i n detail by IWakamatsul 
(|l98lh . Ilnnanen etall (|l983l ) and iBellazzinil (|2004[) . who 
estimate Mq under Newtonian dynamics using the tidal 
radii of remote GCs, the observed tidal radii are affected 
by large uncertainties. The conventional method for deter- 
mining the tidal radii involves a large outward extrapola- 
tion of the surface brightness profile. Innanen et al. (1983) 
suggest random errors of a factor ~ 2 in r^. Moreover, the 
adoption of different models to fit them may lead to sig- 
nificantl y larger limiting radii than e stimated with King's 
models ()McLaughhn fc MevlanI [20031 ). Our first objective 
is to test MOND using the observed tidal radii of remote 
GCs, including fully all the sources of uncertainty in the 
analysis. 

The complete list of known GCs beyond 35 kpc until 
2006 is Pal 2, Pal 15, NGC 7006, Pixis, Pal 14 NGC 2419, 
Eridanus, Pal 3, Pal 4 and AM-1 (|Harrislll996l ). However, 
Pal 2, Pal 15 and Pyxis are not suited for the present 
study because they are strongly affected by extinction. 
Two new extremely low luminosity GCs have been dis- 
covered recently (Koposov et al. 2007), but the number 
of stars detected is not enough to measure precisely their 
luminosities and tidal radii. 

Those GCs with orbital periods, P, much shorter than 
the half-mass internal relaxation time, ir/i, are expected 
to be limited by Galactic tides at perigalacticon, because 
internal relaxation is unable to restore a larger limiting 
radius before the next perigalactic passage (e.g.. Oh et 
al. 1995). In clusters with an orbital period similar or 
larger than the relaxation time, internal relaxation is able 
to repopulate their external regions after perigalacticon 
passage and their tidal radii may be comparable to the 
instantaneous theoretical valu e and, hence, values F ^ 1 
are expected ()Bellazzinill2004D . For the remote GCs of our 
sample from which orbital parameters are available, NGC 
7006 and Pal 3, have P/Uh ~ 1 (|Harris 1996). There- 
fore, we can be certain that they have P ~ 1. Eridanus 
and AM-1, situated at a galactocentric distance of ~ 95 
kpc and 123 kpc, respectively, are likely limited by the 
Galactic tidal force at its present position because their 
trh value ~ 2.5-5 Gyr is short compared to their orbital 
periods. In the case of NGC 2419, trh ~ 35 Gyr. Hence, 
it should be limited by Galactic tides at perigalacticon. 
Finally, P/Uh is difficult to infer for Pal 14 and Pal 4. We 



consider as a working assumption that the observed tidal 
radii of distant GCs, except for NGC 2419, are probes of 
the Galactic tidal force at their present position, bearing 
in mind that we may overestimate rt for Pal 14 and Pal 4. 
The predicted rt values are estimated according to Eq. ([1]) 
for all GCs, except NGC 2419, in a MOND model where 
the Galactic mass is fixed by the inner rotation curve. 

For NGC 2419 we are forced to assume an orbital ec- 
centricity. Since its orbital parameters are unknown we 
will follow a statistical procedure. For the 48 GCs with 
determinatio ns of the orbital par ameters, e = 0.5±0.25at 
the 1(7 level ([Alien. Moreno fc P ichardo 2006). Therefore, 
we will estimate rt for NGC 2419 using Eq. ^ and an 
eccentricity taken as a most likely value of 0.5, with a la 
uncertainty of 0.25. We should note that there is a trend 
on the eccentricity with their apogalacticon distance: most 
GCs with apogalacticon distances < 10 kpc have eccentric- 
ities between and 0.6, whereas those with apogalacticons 
> 10 kpc the eccentricities are populated mainly in the in- 
terval 0.4-1.0. Therefore, adopting e = 0.5 ± 0.25 we are 
overestimating the upper error range on rt since the ec- 
centricity of NGC 2419 is likely > 0.4. 

When estimating rt from Eqs. ([1]) or 0, is derived 
from the distance modulus /iy, the color excess E{B — V) 
and the galactic coordinates, whereas the mass of the clus- 
ter is derived from the apparent integrated 1^-magnitude 
Vt, /iy and the stellar mass-to-light ratio M/L. Redden- 
ing correction, though small (less than 0.1 mag in the color 
excess for all of the GCs considered), was taken into ac- 
count. As said before, in order to quantify the success or 
failure of MOND in explaining the measured tidal radii 
of GCs, we should include the associated uncertainties in 
the analysis: there are dispersions around /xy, Vt^ M/L 
and Keenan's factor. We will proceed adopting the same 
associ ated uncertainties on input parameters as Bellazzinij 
([2004[ ) and are briefly described below. 

(i) The me an values of and Vt were extracted from 
[Harris! ([19960 . The standard deviations for both and 
Vt were taken of 0.1 mag for NGC 2419 and NGC 7006 
and 0.2 mag for the remaining clusters. The errors in /xy 
are thought to account for both the measurement errors 
and the uncertainties affecting the distance scale of GCs 
([Cacciari[[l999l : [Bellazziiii[[2004l ) . No scatter was associated 
to the color excess; it was kept fixed for each GC. 

(ii) The M/L values follow a truncated Gaussian distri- 
bution with a mean value of 1.2 and a standard deviation 
of 0.4, so that M/L = 1.2 ± 0.4. The distribution is trun- 
cated at M/L = 0.5 (e.g., Pryor & Meylan 1993; Feltzing 
et al. 1999; Parmentier fc Gilmore 2001). 

(iii) The MONDian Keenan factor k is extracted from 
a uniform distribution bracketed by the intermediate and 
radial Roche radii, i.e. 0.35 < k < 1.0. In fact, detailed 
analysis of the tidal radii of Galactic GCs in the Newtonian 
casqj suggest that the Newtonian fc values must vary from 
cluster-to-cluster in the range 0.5 < fc < 1.0 (e.g., Bellazz- 
ini 2004), probably because there are several mechanisms 
that induce an outward migration of stars and the escape 
of all extra-tidal stars is by no means assured. Therefore, 
it turns out that the possible k values are bracketed by the 
intermediate and radial Roche radii, which translates into 
0.35 < fc < 1.0 in MOND. 



^ Note that GCs are supposed to contain no dark matter in these studies. 
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The inferred tidal radii with la error bars for Mq = 
0.7 X 10^^ Mq (see ^2.1^ are shown in Fig. lDll For compar- 
ison, the observed Umiting radii are also included. From 
inspection o f the quali t y and extent of the surface bright- 
ness profile. iBellazzinil ( 20041 ) assigns a standard deviation 
for the observed n values of 10% for NGC 2419, 20% for 
NGC 7006 and 30% for the remaining clusters. However, 
estimates of the tidal radii of Pal 3 and AM-1 have been 
superseded by .Hilkerl ()2006f) . Hence we have adopted these 
updated values with an uncertainty of 20%. It can be seen 
that the central predicted values are systematically larger 
than the observed values except for Pal 14. As anticipated, 
however, the uncertainties are so large that it is not possi- 
ble to rule out MOND at all. The most problematic cases 
are Pal 3 and Pal 4, for which there is little overlap be- 
tween the error bars. However, since the observed cutoff 
radii are plagued with systematic errors (e.g., Heggie & 
Ramamani 1995; McLaughlin & Meylan 2003; Bellazzini 
2004), not included in the present analysis, we conclude 
that the present observational data does not allow us to 
distinguish between MOND and dark matter halos. If fu- 
ture more refined measurements simply reduce error bars, 
leaving central inferred values unchanged, tidal radii of 
GCs will indeed rule out MOND. 

Th e diffe rences between our conclusion and those of 
IZhaol ([2005') are probably due to the fact that we included 
all relevant sources of uncertainties in the problem, and to 
having modeled each GC using the observed M/L (with 
their uncertainties), rather than taking a generic value. 

In the following section we shall extend the previous 
analysis of globular clusters to the case of dSph galaxies 
and evaluate the importance of the Galactic tide on their 
structure, which provide a new test for MOND. 

2.2. Dwarf spheroidal galaxies 

dSph galaxies are not necessarily tidally truncated sys- 
tems because the timescale for two-body relaxation, which 
is ultimately responsible for the diffusion of stars to the 
outskirts of the dSph galaxy, may be so large that internal 
relaxation would be unable to smooth any hypothetical 
sharp edge in the stellar luminosity profile in one Hubble 
time. In the case of the remote Leo II dSph, for instance, 
the two-body relaxation timescale is > 300 Gyr even under 
modified dynamics (derived using Eq. 42 in Ciotti & Bin- 
ney 2004 with a ratio between dark matter and luminous 
matter of 7?. « 20) and, hence, the fractional mass-loss 
rate Mj},dMd^/dt < -QMb/Uh ~ -1.5 x lO'^ Gyr"! 
(|Gnedin et al.|[l999l ) , is so small that the photometric con- 
tribution of unbound stars would be gone in the process 
of subtraction of the background population. 

We have selected all the Milky Way satellites placed at 
40 < D < 260 kpc with determinations of their internal 
dynamics: strong candidates of the satellite companions to 
the Milky Way, namely Ursa Major II, Segue 1, Leo T and 
Bootes II dwarves are excluded. Canes Venatici I is also 
excluded because its complex dynamical structure would 
require a more sophisticated analysis (Ibata et al. 2006) 
and Leo IV because the kinematic data available is still 
very poor. The observed limiting radii were taken from 
Irwin & Hatzidimitriou (1995, henceforth IH95). 

The predicted tidal radii under MOND theory were de- 
rived from Eq. ^ using Mq = 0.7 x 10" Mq once Mdw 



is known. We have estimated Mdw in two different man- 
ners. In a MONDian world, the dynamical mass is fully 
determined by the stars. If a value M/L = 3, correspond- 
ing to a normal old stellar population, is assumed for all 
the dSphs, then Fm < 1.25 for all of them except for Sex- 
tans, whose value is Fm = 2.1. This indicates that in 
order to have Fm < 2 for all galaxies, we should invoke a 
M/L ~ 4. Are these values of M/L consistent with their 
internal kinematics? In order to answer this question, the 
MOND M/L values, {M/L)m, were inferred from the in- 
te rnal dyna mi cs of t he stars based on the most recent data. 

iMikroml (|1995D defined a parameter rj = 
1.5{(T/V)^{D/rc), with cr the stellar velocity dispersion 
of the dwarf, Tc its core radius, and V the galactic rota- 
tional velocity at D, which coincides with the asymptotic 
rotation velocity Voo for all the dwarfs. When 77 1, 
the dwarf can be treated as an isolated system and the 
global M/L was obtained exploiting a simple relation be- 
tween the total mass M and the mean line-of-sight velocity 
dispersion, cr, in the isotropic, deep MOND case: 

(jGerhard fc Spergell 119921 : lMilgromllT99l . When 77 < 1 
and all accelerations relevant to the dwarf dynamics are 
smaller than ag, the dynamics becomes quasi- Newtonian. 
The MOND M/L for an object in the quasi-Newtonian 
regime can be obtained immediately from its Newtonian 
estimate, 

'M\ /.gextA ^M' 
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where t/oxt is the acceleration of the system in the exter- 
nal field. We must notice, however, that due to the fact 
that the acceleration becomes anisotropic (Appendix EJ, 
the above relation is not strictly exact. Some idea of the 
uncertainty can be gained from the case of a homogeneous 
sphere, which we derive in Appendix[Bl The mass-to-light 
ratio in that ideal case is underestimated by a fraction 
~ 30% (Appendix |B|) and, hence, due to the uncertainties 
in the structural parameters of these galaxies, this error is 
unimportant. For an isolated spherical system in the deep 
MOND limit, the central, MOND M/L value, (M/L)o,m, 
as a function of the Newtonian estimator {M/L)o,n is 
given by: 

^2 
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with (7n the central velocity dis persion and Tc the core ra- 
dius (iGerhard fc Spergell [199^ . 

In Appendix [C] we briefly describe the observational pa- 
rameters used to infer the mass-to-light ratios for each 
galaxy individually. Bootes, Sculptor, Sextans and Her- 
cules have been assumed to be in the quasi-linear New- 
tonian regime, whereas the rest of the dSph were treated 
in the isolated limit. The derived global {M/L)m and 
{M/L)m are listed, when available, in Table IDll with 
ao = 0.9 X 10"^ cm s"^ (Bottema et al. 2002; see ^TT|) and 
the updated luminosities. For Coma Berenices, Bootes, 
Ursa Major, Hercules and Canes Venaatici II, the val- 
ues reported correspond to the central mass-to-light ra- 
tios. When the derived {M/L)m ratios lie within the 
range 0.5 < {M/L)m < 12, they might be in agreement 
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with a 'naked' old stellar population ijQueloz et al.|[T995t 
iRomanowskv et al] [2003). Only when the resulting value 
is well out this range, we provide a lower limit at the 
2ct confidence level, assuming that dSph galaxies have a 
isotropic velocity distribution; corrections for this are of 
as much as a factor 2. We wish to know if when as- 
suming these lower 2(t limits on {M/L)m, it also holds 
that Fm < 2. However, the reader interested in the mean 
{M/L)m values is referred to Appendix [Ul 

The observed values of the limiting radii inferred by 
IH95 were used to derive and Fm since they are usu- 
ally adopted as the nominal values. The observed rj in a 
certain galaxy may depend on the type of stars selected 
for its determination. For instance, red giant branch stars 
(RGB) are significantly less centrally condensated than the 
blue horizontal branch population in Sculptor (Westfall et 
al. 2006). More importantly, recent studies in some dSphs 
find a "break" in the King profile. The population beyond 
the break radius must signify either the presence of an ex- 
tremely broad distribution of bound stars or the presence 
of unbound tidal debris (e.g., Westfall et al. 2006). In the 
latter case, rt should be identified with the break radius. 
For Ursa Minor, Sculptor, Draco, Fornax and Leo I, we 
have listed additional rows in Table [DT] with new determi- 
nations of rt and luminosity. For conciseness, however, we 
do not give all the possible combinations of rt and L but 
only the most relevant to assess the influence of observa- 
tional uncertainties. 

The derived values of Fm suggest that Leo I is the least 
likely to be affected by tidal effects. On the opposite side, 
Ursa Minor, Sculptor and Sextans may have Fm ^1. Is 
there any evidence of the tidal distortion in these galaxies? 

From photometric data of Ursa Minor, Palma et 
al. (2003) suggest that this object is evolving significantly 
because of the tidal influence of the Milky Way, although 
they cannot determine if the extra-tidal stars are now re- 
ally unbound. In the case of Sculptor, the resultant value 
Fm = 0.8-1.1 is in perfect agreement with its appearance 
(see Section [XT]) . We turn now our attention to Sextans. 
Applying Eq. ^ with A: = 1, M/L = 12, £) = 86 kpc and 
Mg 0.7 X 10" M0, we obtain rt = 2.15 kpc. This value 
should be compared with that reported by IH95 (Table 
ID1[) to obtain Fm = 1.3 ± 0.5, compatible with the re- 
quirement F < 2. We suggest that deep observations of 
Sextans, both photometrically and spectroscopically, out 
to radii larger than those where tails are detectable, may 
provide a test of MOND. If it is confirmed that Sextans 
is not being tidally distorted, this galaxy may need an ex- 
tended halo of dark matter even under MOND dynamics. 
If so, the appeal of MOND would appear to be signifi- 
cantly weakened. In the standard Newtonian dark matter 
scenario -F/v — 1-4 in Sextans, but it is likely that the 
global {M/L)n is larger than the quoted value of ^ 100 
(Kleyna et al. 2004). 

Finally, Draco, which has a large kinematic {M/L)m, 
would survive the Galactic tide even if a {M/L)m — 2, 
typical of a old stellar population, was adopted. 

We conclude that it is very difficult at present to use 
the tidal radii of dSphs to distinguish between MOND 
and standard CDM. Until kinematic measurements defini- 
tively identify the tidal radii, the tidal approach should 
be treated with caution, as it is only an indicator of the 



importance of the Galactic tidal force on the structure. 
In MOND, the dynamical M/L are systematically larger 
than the M/L inferred using the tidal radii (i.e. the re- 
quirement Fm < 2). Therefore, explaining the internal 
dynamics is a more profound question than the tidal radii. 
In fact, the high {M/L)m ^ 20 values required in Will- 
man 1, Coma Berenice, Ursa Minor, Draco, Ursa Major 
and possibly Bootes are difficult to accept if MOND is a 
valid alternative to dark matter in dSph galaxies. Will- 
man 1 , Coma Berenice and Ursa Major are unique targets 
to test the success or failure of modified gravities at galac- 
tic scales. One could argue that an appropriate choice of 
anisotropy in the orbits of stars within these galaxie s could 
reduce somehow these estimates. However. iLokasI (|2002f) 
found the opposite trend in Draco and Fornax. 

The relation between M/L and luminosity is shown in 
Figure [D2l for our sample of galaxies. This figure is based 
on that given in Mateo et al. (1998) but for MOND. Even 
with considerable scatter, the clear trend between M/L 
and luminosity is because MOND is based on a charac- 
teristic acceleration and dSph galaxies present a charac- 
teristic velocity dispersion. According to this relation, the 
kinematics of the faintest dSph galaxies will be puzzling 
under MOND. 

Galactic tides could, in principle, affect the kinemat- 
ics and the apparent M/L derived from kinematic stud- 
ies that assume dyn amic al equili b rium. For instance, 
iKuhn fc MiHal (|1989[ ) and lKroupal (|1997D presented tidal 
models to explain large values of M/L in dwarf galaxies 
without restoring to dark matter. However, there is grow- 
ing evidence that galactic tides cannot infiate the global 
M/L values: tides produce large ordered motion s rather 
than large random motions (jPiatek fc Prvoij|l995f ). Oh et 
al. (1995) found that the unbound but not yet dispersed 
systems have velocity dispersions that are comparable to 
the virial equilibrium value prior to disruption. Therefore, 
if the measurements do not seriously overestimate the ve- 
locity dispersion, either the mentioned galaxies are not 
in dynamical equilibrium due to recent gas mass loss or 
to a major merger with another dwarf galaxy, undergoing 
still a relaxation process, or they contain large amounts 
of dark matter even under modified dynamics. We believe 
that it is unlikely that the gas has been stripped from 
these galaxies only very recently, in the last few internal 
crossing-times. In particular, the fact that Draco and Ursa 
Minor have only very old stellar populations (> 10 Gyr), 
suggests that they have possessed a dynamically-negligible 
fraction of mass in gas since a long time ago (Mayer et 
al. 2007). 

3. RETENTION OF LOW-MASS X-RAY BINARIES 

Recently, iDehnen fc Kind (|2006f ) have suggested that 
in order for Sculptor to retain the population of low- 
mass X-ray binaries observed by Maccarone et al. (2005) 
an extended dark matter halo of > 10^ Mq is required, 
i.e. {M/L)m ~ 600. Here we explore the implications of 
this observation in order to see if the stronger gravity due 
to the stars alone under MOND also has the ability of 
retaining the observed LMXBs or not. 

3.1. The Sculptor dwarf spheroidal: Internal parameters 
and proper motion 
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The Sculptor dSph is a satellite galaxy of the Milky Way, 
which presents a luminosity (1.4 ± 0.6) x 10^ L© (IH95). 
According to the proper motion measurements by Piatek 
et al. (2006), it is on a polar orbit with apogalacticon at 
122 kpc and perigalacticon at 68 kpc. The current galac - 
tocentric distance of Sculptor is 79 ± 4 kpc (|Mateolll998D . 
A King model profile of limiting radius 79'. 6 ± 3'. 3, cor- 
responding to 1.85 kpc and core radius 7'. 14 ± 0.'33 (165 
pc), fits the density profile well out to 60'. Beyond this 
radius, a break population, which extends up to 150', has 
been found (Westfall et al. 2006). The ellipticity of the iso- 
density contours increases with increasing projected radius 
from the center of the dSph, from rounder to a value 0.32 
in the outermost region (IH95; Westfall et al. 2006). Since 
there is no evidence for rotation, this flattening should be 
due to anisotropy in the velocity distribution, but it may 
also indicate that Sculptor is tidally distorted. 

There is evidence that Sculptor contains two distinct an- 
cient (both > 10 Gyr old) stellar components. Tolstoy et 
al. (2004) derive different velocity dispersions for Sculptor 
stars separated into metallicity groups. The line-of-sight 
velocity dispersion of the metal-rich component is ~ 7 ± 1 
km s~^ , which is similar to the central line-of-sight velocity 
dispersion in Sculptor 6.2±1.1 km s~^ (Queloz et al. 1995; 
Westfall et al. 2006). Interestingly, the metal-poor com- 
ponent has a velocity dispersion of w 11 ± 1 km s~^ (see 
also Clementini et al. 2005). Westfall et al. (2006) derive 
a global velocity dispersion of 8.8 ± 0.6 km s~^, which 
is intermediary to those that Tolstoy et al. (2004) report 
for these mctallicity-separated populations. In the follow- 
ing, we will take the latter value as the mean line-of-sight 
velocity dispersion within the core, Uc- From Fig. 12 in 
Westfall et al. (2006), one can see that it is a generously 
taken value. 

3.2. Escape velocity in the quasi- Newtonian MOND 
limit: Sculptor as a reference case 

The effective gravitational mass of a isolated galaxy un- 
der MOND increases linearly with i?, and hence, any test 
particle is bound to the potential. However, due to the 
nonlinearity of the MOND field equation, for a satellite 
system situated in an external gravitational field, there ex- 
ists a radius at which particles are stripped from the dSph 
by the tidal field of the host galaxy (see For certain 
binary kick velocities, the LMXB system could become 
unbound from the gravitational potential of the satellite. 
Our aim in this Section is to estimate the kick velocity of 
a presupernova binary necessary to acquire the escape ve- 
locity, denoted by vu- In Sculptor, the calculation greatly 
simplifies because of the following reasons. It can be seen 
that the internal acceleration ~ 1.6cr^/rc ~ 2 x 10~^ cm 
s~^ felt by a star in the core of Sculptor is significantly 
smaller than aq and thus, for the dynamics of an object 
in Sculptor the deep MOND regime would apply. In addi- 
tion, the parameter »7 is < 0.8 for V > 185 km s^^. Since 
the internal acceleration goes as oc r at small radii and as 
< r^^ at large radii, this value of 77 points to domination of 
the external field throughout the dwarf galaxy for V > 185 
km s""'^ (see also Milgrom 1995). At V = 170 km s~^, we 
get ?7 ~ 1 , implying it is in the borderline of isolation with 
respect to the external acceleration field. Neither limit is 
valid in this case but the two limits should give similar 



values for the escape speed. 

A rough derivation of the escape speed follows; a more 
detailed calculation is presented in Appendix [Dl Suppose 
that a certain dSph can be described by a one-component 
King model (i.e. mass follows light). Then, the mass pro- 
file can be parameterized by Wq, po and Tc, with po the 
total mass density at the center. In Newtonian dynam- 
ics, Wq is the depth of the potential in units of the square 
velocity-dispersion parameter ct^ = AnGporl/9. The es- 
cape velocity at any position is, therefore, smaller or equal 
to the central escape velocity \/2Woa^. For the Galactic 
dSphs, Wo lies in the range 2 < Wq < 5. 

Consider now what happens under MOND theory. If 
this galaxy is in the quasi-Newtonian regime and MOND 
has the ability of reproducing the dynamics of this galaxy 
without dark matter, the depth of the potential under 
MOND must be similar to that in the case of the Newto- 
nian dark matter scenario. Hence the central escape veloc- 
ity is also \/2Wo(j^ under MOND. In the case of Sculptor 
W^o ~ 4 (IH95; Walcher et al. 2003; Westfall et al. 2006). 
For this value of Wo, the velocity-dispersion parameter is 
related to the observed velocity dispersion within the core 
through (7 = 1.16crc (see Fig. 4-11 in Binney & Tremaine 
1987). Adopting CTc = 8.8 km s~^ (see t|3.ip . the central es- 
cape velocity in Sculptor under MOND is ~ ^/l(^8a^ = 29 
km s~^. 

The escape speed depends, of course, on radius. For 
instance, the escape velocity at the core radius of the 
same King model has been calculated to be 23 km s~^. 
Since the rms velocities of the progenitors in the cen- 
tral parts of Sculptor are ^ 8.8 km s^^, kick velocities 
of V232 - 8.82 _ 21.5 km s'^ are enough for a LMXB 
formed at a core radius to be expelled. Consequently, if 
the centers of mass of all LMXBs have a recoil velocity 
of 21.5 km s~^, a fraction of them are retained and the 
rest becomes unbound, depending on the location where 
the kick occurs. Since the fraction of mass within the core 
radius is 35%, a significant fraction (~ 65%) of the new 
formed LMXBs will escape for a kick velocity of 21.5 km 
s~^ (a statistical approach to infer the mean kick velocity 
is given in Appendix iDl reaching the same result). 

The systems that survive as binaries and become LMXB 
progenitors attain a system velocity of 180 ± 80 km s~^, 
if the distribution of angles between the kick velocity and 
the orbital plane of the presupernova binary is isotropic. 
If the range of kick directions is restricted to a cone along 
the spin axis with an opening of 20° , LMXBs are launched 
with kick velocities of order 20-100 km s~^ (e.g., Podsi- 
adlowski et al. 2005). In this case and if the distribution 
of kick velocities were drawn randomly in the range 20- 
100 km s~^, less than 2 percent of the LMXBs would be 
retained in Sculptor. The recent analysis of the Galactic 
pulsar proper-motion data by Arzoumanian et al. (2002) 
found evidence for a bimodal distribution, with ^ 40% of 
the pulsars contained in a Maxwellian component with a 
dispersion of ~ 90 km s~^, and the remaining neutron stars 
in a Maxwellian component with a dispersion of ^ 500 km 
s^^. With this probability distribution, the retention frac- 
tion is slightly larger but still less than 2.8%. 

The presence of a heavy binary companion at the time 
of the supernova explosion will make the retention more 
hkely (e.g., Davies & Hansen 1998). Pfahl et al. (2002) 
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consider the inclusion of binaries and conclude that the 
retention fraction is probably not larger than several per- 
cent when they apply a single Maxwellian fast kick mode 
at 200 km s~^ and a central escape speed of 50 km s~^. 
Scaling down their results, this implies that the retention 
fraction in Sculptor is likely not larger than a few per- 
cent even assuming that the velocity kick distribution is 
described by a slow Maxwellian mode at about 90 km s~^. 
It is therefore puzzling to understand how Sculptor was 
able to hold on to most of its X-ray binaries in MOND 
theory. 

One possibility is to invoke dark matter in MOND. Even 
if Sculptor obeys MOND, it can also contain a dark mat- 
ter halo more extended that the optical galaxy such that 
the kinematics within the visible extent of Sculptor dwarf 
would not be greatly affected, but will be very effective at 
holding on to the LMXBs in the dSph potential. In this 
model, the MOND paradigm lowers the discrepancy be- 
tween the binding mass and the baryonic mass but it still 
requires a total mass-to-light ratio (M/L ) m 120, where 
we ha ve used Eq. @ with (Af /L) n > 600 (i Dehnen fc Kind 
[200^ . A gain, the inclusion of a dark halo in MOND weak- 
ens its appeal. 

3.3. Compact early conditions 

It is simple to see that if the stellar distribution was sig- 
nificantly more compact in the past than it is at present, 
the probability of old LMXBs to be bound to the poten- 
tial increases. This is true for both Newtonian and MOND 
dynamics, but applies only if the compact dSph was subse- 
quently heated to its current condition by a slow adiabatic 
process. A dSph galaxy could have been more massive in 
the past since mass-loss events can be caused by ram pres- 
sure stripping, supernovae explosions and tidal stripping. 
However, since these mechanisms are thought to be sud- 
den and violent, it is very unlike ly that they are a ble to 
produce the desired effect. While iRead fc Gilmord (|2005f ) 
have simulated the dynamical effects of mass loss on the 
remaining stars and dark matter in the Newtonian case, 
it is illustrative to see that, in fact, the velocity disper- 
sion of baryons is not altered enough to consider mass loss 
as a promising alternative possibility to explain the re- 
tention of observed LMXBs. In the recent simulations by 
Mayer et al. (2007), the formation of the darkest dSphs 
like Draco and Ursa Minor, is the result of the transfor- 
mation from a gas-rich dwarf to a dSph by repeated tidal 
shocks at pericenters when the gas is readily removed by 
ram pressure stripping. Again, these impulsive shocks are 
not able to change the stellar properties in an adiabatic 
way. In this scenario, dSphs that fall into the Milky Way 
halo late, suffer little central stripping and may have peri- 
odic bursts of star formation because of tunneling of this 
gas by a tidally-induced bar. 

In the dark matter picture, there might exist another 
somewhat speculative via to produce dynamical heating 
of the stellar population. If the dark matter component 
is made up by massive black holes, gravitational encoun- 
ters with stars will produce a transfer of energy between 
black holes and stars, in the attempt to reach equiparti- 
tion. iJin. Ostriker fc Wilkinsoiil ()2005f ) have demonstrated 
that the dynamical heating of an initially compact stellar 
distribution might produce a remnant distribution similar 



to those in dSphs for black hole masses in excess of 10^ 
Mq. LMXBs formed at early stages (> 3 Gyr ago) could 
become bound to the potential well of the dSph galaxy. 
However, this scenario is highly unlikely (Hernandez et 
al. 2004; Sanchez-Salcedo & Lora 2007). 

3.4. Comparison with the escape velocities within the 
dark matter paradigm 

For a central line-of-sight velocity dispersion of 6.2± 1.1 
km s~^ as measured in Sculptor, the central mass-to- light 
ratio in the classical Newtonian view was estimated to be 
6-13 (Queloz et al. 1995; WestfaU et al. 2006). Given its 
large uncertainties, this range of values is consistent or 
marginally larger than the stellar mass-to-light ratio ex- 
pected for an old stellar population (Queloz et al. 1995). 
Nevertheless, its kinematics could be also compatible with 
Sculptor being hosted by an extended subhalo. In fact, 
a mass estimator for dark matt er halos via the escape - 
velocity argument was used by iDehnen fc Kind ()2006[ ). 
They construct diagnostic dark halo models with cumula- 
tive mass profiles described by a simple analytic formula, 
designed to explore the generic parameters required for 
retaining LMXBs in Sculptor. These halos have different 
core radii, but all satisfying that the dark mass within the 
visible galaxy, with a radius of ~ 1.85 kpc, be 5 x 10^ M©. 
They find that no reasonable dark matter halo in Sculptor 
can retain all LMXBs with ejection velocities in the range 
20-100 km s'^. To hold on to LMXBs with ejection ve- 
locities of up to 60 km s^^, a dark matter halo having a 
total radius of upwards of 15 kpc is required. LMXBs with 
lower ejection speeds are more easily retained and require 
dark matter halos with more conservative parameters. 

In order to make the comparison with MOND as fair 
as possible, we have explored the parameter space that 
ensures that LMXBs are retained in Sculptor but using 
full King models for the dark halo of Sculptor . Consis- 
tent with the findings of iDehnen fc Kind ()2006[ ). we find 
that no realistic halo model is sufficient to retain LMXBs 
having large ejection velocities. Given that these ejection 
velocities are expected in the range of 20-100 km s~^ (e.g., 
Podsiadlowski et al. 2005), and that the velocity disper- 
sion in Sculptor is of only around 9 km s^^, this conclu- 
sion is not surprising. Actually, for the dark mass limits 
within the visible galaxy coming from dynamical studies, 
and imposing halo total radii smaller than 15 kpc, little 
room for varying the King halo parameters remains. For 
reasonable parameters of the dark halo, escape velocities 
are always of order 30-40 km s~^. Figure lD3l shows one 
such mass model for Sculptor, with the upper panel giv- 
ing the cumulative mass profile for the stars (also taken 
as a King distribution with a stellar mass-to-light of 1.5 
and a concentration of 0.63), and the total mass. The 
dominance of the dark matter halo is dramatic, especially 
at large radii, although not extraordinary in comparison 
to inferences of large galaxies. The bottom panel gives 
the escape speed as a function of radius together with the 
components due to the stars and dark matter. We see that 
even at the very center, the escape velocity is of only 40 
km s~^, and drops to 30 km at a distance of around 
0.3 kpc. One should perhaps expect the progenitors of 
LMXBs to have been located within the central 0.165 kpc 
of Sculptor, which defines the core radius of the stellar 
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populations. This leads to the conclusion that standard 
dynamics and a reasonable dark matter halo proposal for 
Sculptor can explain the retention of LMXBs, only for 
those having started with a relatively low < 40 km 
initial kick. This is a more satisfactory situation than in 
the MOND case. In order to keep neutron stars with kick 
velocities > 50 km bound, an extraordinarily extended 
and massive halo with a very large core radius should be 
invoked. Additionally, if a large fraction of neutron stars 
in dwarfs were formed in binary systems, the retention of 
LMXBs in Sculptor could be more plausible, but not triv- 
ial, certainly not for all initial velocities in the theoretical 
ejection range of velocities, at least not with the restric- 
tion that the total M/L < 35 within 1.85 kpc. Therefore, 
the problem of explaining the large retention of LMXBs is 
not exclusive of GCs (e.g., Pfahl ct al. 2002). 

4. CONCLUSIONS 

We have pushed the measurements of tidal radii in 
Galactic GCs and dSph galaxies attempting to distinguish 
between dark matter and MOND. Except for Pal 14, the 
predicted tidal radii of GCs are systematically larger than 
the observed nominal values. However, this is not enough 
to rule out MOND because after including properly all 
the uncertainties, we find that they are consistent with 
the observed values at a Icr level. The importance of the 
Galactic tide on the survival of dSph galaxies has been 
also investigated. Assuming mass-to-light ratios compat- 
ible with a naked stellar population, the present dSph 
galaxies preserve their integrity except perhaps Sextans, 
which might be undergoing tidal disruption. Wc suggest 
that deep observations of Sextans might provide a further 
test of MOND. 

Based on the most recent data and with the updated 
values of ao and Mq, the dynamical mass-to-light ratios 
inferred under Newtonian and MOND gravities have been 
derived and compiled to sec if new data can definitively tell 
us what the law of gravity at small galactic scales is (~ 300 
pc). Since the dynamical mass in MOND should be domi- 
nated by stars, it is worrying that it requires mass-to-light 
ratios > 12 for those dSph galaxies with radially extended 
data right to the optical edge (Ursa Minor and Draco). 
We warn that Willman 1, Coma Berenice, Ursa Major and 
Bootes may be also problematic for MOND. In particular, 
preliminary estimates for Willman 1, Coma Berenice and 
Ursa Major suggest M/L > 20, at the 95% confidence 
level, even when possible outliers that inflate the inferred 
velocity dispersion, and hence the mass-to-light ratio, are 
omitted. That MOND requires a fraction of dark matter 



in clusters of galaxies is a well-established issue. Taken 
at face value, our results indicate that MOND needs also 
a dark component in dSph galaxies with a mass fraction 
much beyond that required in galaxy clusters, weakening 
the appeal of the MOND paradigm, although this fact by 
itself does not rule out MOND as the gravity law at low 
accelerations. 

An estimate of the mass of Sculptor has been given via 
the escape- velocity argument: LMXBs are ideal probes 
of the total mass as they should penetrate well outside 
the visible galaxy due to their high recoil velocities > 20 
km s~^. In the standard dark matter scenario, LMXBs 
with velocities > 40 km are difficult to be retained 
for reasonable parameters of the dark halo. Hence, the 
well-known retention problem of LMXBs in GCs persists 
in dSphs. In the MONDian case, we have calculated that 
the mean square kick velocity of a prcsupcrnova binary 
to acquire the escape velocity is ^ 24 km s^^ in Sculp- 
tor. About 95% of the formed LMXBs would have escaped 
and been stripped from the Sculptor galaxy. In order for 
Sculptor to hold on to LMXBs with kick velocities of 50 
km s~^, a mass-to-light ratio > 120 is required even un- 
der MOND hypothesis, making it structurally different. 
Therefore, either this is new evidence for the necessity of 
an extended dark halo in a galaxy in MOND, or the con- 
ventional thinking regarding neutron stars kicks must be 
modified. Measuring the radial velocities of the observed 
LMXBs would provide an obvious test of alternative grav- 
ities. High radial velocities > 30 km s~^ may confirm the 
presence of a very massive dark matter halo much more 
extended than the stellar population. 

Finally, one could argue that the large values of the 
M/L in MOND must be attributed to a problem with the 
velocity dispersion measurements and their estimated un- 
certainties. The required change in the velocity dispersions 
appears unlikely on purely statistical grounds and, in ad- 
dition, they do not act in a systematic manner in Draco 
and Ursa Minor when compared to Sculptor. 

In a MOND universe, data for several of these systems 
imply the presence of an extended dark matter halo. Al- 
though not a disproof of MOND, it does reduce the appeal 
of MOND significantly. 

We would like to thank the anonymous referee for a 
carc;ful reading of the manuscript and useful suggestions 
which improved the presentation of our final version. The 
authors acknowledge support from PAPIIT project IN- 
114107-3. 



APPENDIX 

MOND IN THE QUASI-NEWTONIAN LIMIT: VIRIAL THEOREM 

For a satellite galaxy embedded in a constant external field ^ext) the gravitational field equation that governs the 
kinematics of the stellar component in the quasi-Newtonian limit is 

V2 + $ = ^4nGp. (Al) 

The solution of this modified Poisson equation for a distribution of mass p{-r) is given by 



= _ ( ^ W ^^r- (A2) 

' ffext y J ^/2(x' - x)2 + 2(y' - J/)2 + (z' - zf 
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As usual, the virial theorem for a stationary system can be derived after integration over an arbitrary volume V: 

pv'^(fif= / pr-^^(fr. (A3) 
V Jv 

Substituting Eq. into Eq. ((X3|) . one gets: 

2K + W^Q, (A4) 
where K is the total kinetic energy of the system and W is the system's total potential energy: 

J p(f)$(f) (ff. (A5) 

We have seen that the virial theorem has the same form than in the Newtonian case, with $ as given by Eq. (|A2p . 

POTENTIAL ENERGY OF A HOMOGENEOUS SPHERE IN THE QUASI-LINEAR LIMIT 

We wish to calculate the potential energy, 14^, of a homogeneous sphere of constant density p and radius vq in the 
quasi-linear MONDian limit. According to Eq. (jASp . we need $ in order to evaluate W. The field equation (|A1[) . can be 
transformed into the standard Poisson equation by making the substitution z' = z/^/2: 

V2$' = ^4^Gp(x, y, V2z'), (Bl) 

5cxt 

where $' = ^'{x,y,z') and ~ d^/dx^ + d^/dy^ + d^/dz'^. The relationship between the potentials is ^{x,y,z) = 
$'(a;, y, z/\/2). For a sphere of constant density, satisfies: 

V2$' = ^47rGp', (B2) 

5oxt 

where p' — p within the ellipsoidal body bounded by the surface Tq = i?^ + 2z'^, which is an oblate ellipsoid with 
eccentricity e = 0.71. We are now in a position to calculate W — 1/2 J p'^'cfif' taking advantage of the theory of 
homoeoids (e.g., Roberts 1962; Binney & Tremaine 1987): 

with /(e) an analytic function of the eccentricity as given in Table 2-2 of iBinnev fc Tremaind (|1987D . being J = 1.57 for 
e = 0.71. 

Finally, we can obtain W, W — 1/2 / p^(Pf, in terms of W' by expressing W in terms of the new variable of integration 
z ^ z' . This leads to = \/2W' . Substituting the value of W' from Eq. (|B3p as a function of the total mass M — AnprQ/i, 
we find 

W = -0A7(^]^. (B4) 

If the MOND dilation along the {x, y)-directions were ignored and the potential were assumed to be isotropic rather than 
anisotropic, the potential energy would be — 0.6(ao/(?ext)GM^/ro. Hence, the error made using the isotropic approximation 
is 28%, for a homogeneous sphere. 

DYNAMICAL MASS-TO-LIGHT RATIOS OF DSPHS 

The closest and faintest dwarf in our sample is Willman 1 at a distance of 40 kpc and with a luminosity of ~ 855 Lq 
(Martin et al. 2007). For the reported core radius of 16.6 pc, this galaxy may be treated as isolated for internal velocity 
dispersions > 2.7 km s^^, corresponding to masses of 8.5 x 10"* Mq. Therefore, for realistic mass-to-light ratios of a stellar 
population, Willman 1 is in the quasi-linear regime. The expected velocity dispersion is: 




OMJ^J—, (CI) 



where r^p is the projected half-mass radius (e.g., Baumgardt et al. 2005). Suppose for a moment that the mass-to- light 
ratio is ^ 3. Under this assumption, Eq. (jCip implies a = 0.53 ± 0.1 km s~^, where rhi, « 21 ± 7 pc was used. By 
contrast, Martin et al. (2007) measure a radial velocity dispersion of 4.3±^-3 km s~^ and higher than 2.1 km at the 
99% confidence limit. Therefore, if Willman 1 does not contain dark matter, it is already gravitationally unbound even 
under MOND, and will disperse in a few crossing-times, i.e. in a few ~ 10'' yr. This is at odds with the photometric data 
that do not show tidal tails or visible increase of the velocity dispersion with distance (Martin et al. 2007). Thus it would 
seem that Willman 1 needs a dark matter halo even in MOND. 

In order for Willman 1 to be in virial equilibrium with a velocity dispersion of ~ 4.3 km s^^, we require a mass-to- light 
ratio between 370 and 550. This estimate was derived from the Newtonian values reported in Martin et al. (2007) but 
noting that since the system is not in deep MOND in this case (the internal accelerations a re of the order of aq), the 
MOND M/L value is only a factor of ~ 1.3 smaller than the Newtonian values (Milgrornl ll994[ ). Adopting the lower limit 
of 2.1 km at the 99% confidence level, Willman 1 lies in the quasi-linear deep-MOND regime and then {M/L)m > 30. 



10 



For the new Galactic dSph candidate found in the consteUation of Bootes, Martin et al. (2007) have derived a velocity 
dispersion B.Sif-g km s~^ from 24 members of the Bootes dwarf. For the adopted distance of 60 kpc, this galaxy has 
a core radius of ~ 225 pc and is well in the quasi-Newtonian regime. Adopting a total luminosity of 1.6 x 10''' Lq to 
8.6 X lO'* L0, implies (M/L)o,Af ~ 120 to 23. Given the uncertainties in the velocity dispersion, (A//i)o,M ^ 15 at the 
68% confidence. We must stress here that the magnitude and associated errors in the velocity dispersion and structural 
parameters of Bootes and Willman 1 are currently crude estimates. 

For Ursa Minor, Wilkinson et al. (2004) derived a line-of-sight rms dispersion of 12 km s^^, with a > 9.5 km s^^ 
at the 95% confidence level. A robust estimate of the total mass can be found from Eq. ([3]) provided that the velocity 
distribution is isotropic. For a value of cr = 12 km s~^ we get {M/L)m = 60 even taking the luminosity derived by Palma 
et al. (2003), while {M/L)m > 23 at the 95% confidence level. 

Based on the Newtonian value {M/L)n = 9 and in the quasi-Newtonian limit (Eq. [5]) with gcxt — V'^/D and V = 170 
km s~^, the global {M/L)m in Sculptor is 1.0 ± 0.5, for the observed central line-of-sight velocity dispersion 6.2 ± 1.1 
km s~^ (Queloz et al. 1995; Westfall et al. 2006). The real global M/L ratios are probably a factor x2 larger because 
the velocity dispersion shows a rise, leading to a mean velocity dispersion in the core of 8.8 ± 0.6 km s""' (Westfall et 
al. 2006). In any case, the mass-to-light ratios in MOND are perfectly consistent with an old stellar population. In a 
Newtonian dark matter view, a flat velocity-dispersion profile well beyond the several core radius, as occurs in Sculptor, 
is interpreted as the presence of a halo that begins to dominate dynamics increasingly outwards. However, in the MOND 
quasi-Newtonian limit, such a behavior for the velocity dispersion, even if the actual values of the stellar mass-to-light 
ratios remain within plausible ranges, would imply a radial change in the intrinsic properties of the stars, something hard 
to explain in such small systems, or t hat the dwarf galaxy is being heated by Galactic tides. 

Draco and Fornax were studied bv iLokai (|2001l [20021 ). The {M/L)m values for isotropic models after updating by 
flo and the luminosities are 29 ± 9 and 2.2 ± 0.4, respectively. The mass of Leo II and Leo I have been estimated with 
Eq. ([3]) using the mean velocity dispersion 6.6 ±0.7 km s~^ (Koch et al. 2007b) and 9.9 ± 1.5 km s~^ (Koch et al. 2007a), 
respectively. For Carina, we estimate {M/L)m with Eq. ^ adopting the central value 6.97 ± 0.65 km (Munoz et 
al. 2006a). 

For Sextans, the two-component models by Kleyna et al. (2004) suggest {M/L)m between 60-300 within 1 kpc, being 
the exact value sensitive to the form assumed for the outer fall-off in the velocity dispersion profile. Based on these 
Newtonian values, assuming that the external field is dominant in Sextans, a {M/L)m of ^ 7.5-36 was found when 
adopting a distance of 86 kpc. Although M/L of GCs have a 3cr range of 1-4.3 (Feltzing et al. 1999; Parmcntier & 
Gilmore 2001), values as high as 7-8 cannot be ruled out for a stellar old population (Queloz et al. 1995; Romanowsky 
et al. 2003). Hence, a mass-to-light ratio of ~ 8 seems to be compatible with its kinematics and might be associated to a 
naked stellar old component. 

In the case of the recently discovered Ursa Major dSph, for a central velocity dispersion of 7.6 km s~^ and a core radius 
Tc ~ 200 pc (Simon & Geha 2007), we obtain a central mass-to-light ratio of 135 in the isolated limit. For the value 
reported by Martin et al. (2007), a = 11.9±2 | km s~^, the mass-to-light ratio is even larger. The velocity dispersion is in 
fact greater than 6.0 km s~^ at a 95% confidence level. Adopting this lower value, Ursa Major lies in the quasi-Newtonian 
regime and, then, (M/L)o,m > 52 with 95% confidence. The magnitude of Ursa Major is currently a crude estimate. 
Improving the data on this galaxy is key to exploring how MOND behaves at the smallest galactic scales (~ 300 pc). 

Measurements of the velocity dispersions for Coma Berenices, Hercules and Canes Venatici II were reported in lSimon h Gehal 
([2007). They measure the velocities of 59 stars in Coma Berenice, obtaining a velocity dispersion 4.6 ±0.8 km s ^ . Equa- 
tion ([5]) with Tc — 41 pc and {M/L)^ — 448 implies a MOND mass-to-light ratio of 105. Hercules is likely dominated by 
the external acceleration because rj ~ 0.91 for Tc = 205 pc and assuming a distance of 140 kpc. For a velocity dispersion 
5.1 ± 0.9 km s^^ and {M/L)m — 332, we obtain {M/L)o^m = 26±9^. Finally, the ultra-faint Canes Venatici II is located 
at 150 kpc. We have used Tc = 84 pc and a = 4.6 ± 1 km s~^ to infer (Af/L)o.M = 38.5±24. 

Since the global luminosity is less certain than the central luminosity, we can use the central mass-to-light estimators 
(M/L)o,M as a check of the robustness of the estimates, for Draco and Ursa Minor. Assuming that mass follows light, 
Odenkirchen et al. (2001) obtain (M/L)jv — 146 for Draco. If we take this value as representative of the Newtonian value 
within the optical radius of Draco, and for a central velocity dispersion ctq = 8 km s~^ (Lokas et al. 2006) and core radius 
of 180 pc, Eq. (P gi ves (M / L)n,Ar — 240. In the case of Ursa Minor, the central Newtonian value is 100 after updating 
the quoted value of iMateol fl998l w ith the latest determination of ctq = 12 km s~^ (Wilkinson et al. 2004), instead of 
9.3 km s^^ adopted in lMateol |l998D . From Eq. ([5]) with Tc = 200 pc, we get (M/L)o,m = 33. These central values of 
(M/L)o,M are therefore in agreement with the global estimates. 



MEAN SQUARE KICK VELOCITY TO ESCAPE 

Consider a progenitor system of a LMXB, which is embedded in a galaxy that lies in the quasi-Newtonian regime in 
MOND. If the progenitor system is moving at velocity v and the supernova explosion occurs at r, then the minimum kick 
velocity for the new formed LMXB to become gravitationally unbound, Vk, satisfies: 

i [v + vkf + $(f) -^\nxr\^ = $(fO - ^1^^ X nf, (Dl) 



^ IMilgromI l|1995f) obtained a central M/L = 4, using ctq = 9.2 km s ^ and ao = 2 X 10 * cm s ^. 
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where fl is the angular velocity of the dwarf satellite around its parent galaxy, ft = (0, 0, rt) and rt the tidal radius along 
the line pointing the centers of the galaxies. The mean square kick velocity for stars to escape is: 

(vl) = - (v^) - 2 mf)) - ^n j) + < |fl X rf > -\n X rip, (D2) 
where the brackets represent the mass average: 

<«>^^^^^/*>^'^^ 

In Eq. (|D2p we used {v ■ Vk) = for an isotropic distribution of kick directions. 

In the quasi-Newtonian MOND regime. Equation (|D2p can be expressed in terms of observable quantities as follows. 

First, the quantity (\Q x is exactly (r^) /3 for a spherical distribution of matter. Corrections for non-sphcricity 

of the distribution in the central parts of the dwarf satellite were found to be unimportant as far as (f ^) concerns. Ignoring 
the variation of the external gravitational field in the vicinity of the satellite, {v'^'^ can be written: 



«> ^ -liJ^-s;;; U<'->-^'j' 



(D4) 



where K is the total kinetic energy of the stellar system and W the total potential energj0 (Appendix |^ . It is easy to 
show that, independently of the concentration, ^r^) < r| /5 in any King model. Using the virial theorem in the quasi-linear 
regime (see Eg. IA4I) . with K — SMdwO''^ /2, a the global one-dimensional velocity dispersion, fl — V/D and approximating 
$(n) ~ -~G{ao/gcxt)Mdw/rt, it holds that: 

The second term of the RHS can be expressed in terms of a exploiting again the virial the orem and n oting that the potential 
energy in a King model with Wq = 4 is \W\ — (9/5)GM?/rt in Newtonian dynamics (|Kinej|l966l ). which translates into 
an upper limit in MOND due to the dilation factor fEa. lA2p . i.e. \W\ < {9/5)G{ao/gcxt)M^^/rt. Consequently, 

(-^)<8a^-|(^)V^. (D6) 

For the parameters of Sculptor: cr ~ O.OSctc = 8.4 km s^^, rt = 1.8 kpc, V — 170-185 km s""'^, D ~ 79 kpc, we obtain 
(^2)'/' < 23.5 km s-i. 
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Table D1 

A SUMMARY OF RELEVANT PARAMETERS FOR THE CLOSEST dSpH GALAXIES. THE SECOND (AND THIRD) 
ROW FOR SOME OF THE OBJECTS GIVES THE VALUES REPORTED BY OTHER AUTHORS. 



Galaxy 


L 

lO^Lo 


rt 

arcmin 


(M/L)^r- 


Comp." 


Fn 


(M/L)m'= 


Fm 


References for 
(L,n,(M/L)jv) 


Willman 1 


0.0085 


~ 7.4 


500 


1 


0.25 


> 30(95%) 


0.36 


M07, M07, M07 


Coma Ber. 


0.027 




448±ii 


1 




> 46(95%) 




B06, — , S07 


Bootes 


0.86 




130 


1 




23±J^ 




M06b,— , M07 


Ursa Minor 


2.9 


50.6 ±3.6 


79 ±39 


1 


0.55 


12 ±6 


0.51 


M98, IH95, M98 




5.8 


77.9 ± 8.9 


40 ±20 


1 


0.85 


6±3 


0.77 


P03, P03, P03 




5.8 


~ 34 


340 ± 240 


2 


0.19 


> 23(95%) 


0.16 


P03, W04, W04 


Sculptor 


14 


79.6 ± 3.3 


9±6 


1 


1.13 


1.0 ±0.5 


1.10 


IH95, Wc06, Q95 




14 


- 60 


9±6 


1 


0.86 


1.0 ±0.5 


0.83 


IH95, We06, Q95 


Draco 


2.6 


28.3 ± 2.4 


> 240 


2 


0.24 


> 13(95%) 


0.28 


M98, IH95, KOI 




2.4 


> 40.1 


140 ± 40 


1 


0.42 


> 12(95%) 


0.40 


OOl, OOl, OOl 


Sextans 


5.0 


160 ± 50 


107 ± 72 


1 


1.40 


22 ± 15 


1.05 


M98, IH95, IH95 


Carina 


4.3 


28.8 ±3.6 


43±fi 


1 


0.42 


4.5 ±0.9 


0.35 


M98, IH95, M06a 


Ursa Major 


0.15 




1024 ± 636 


1 




> 52(95%) 




BOO — , S07 


Fornax 


155 


71.1 ±4.0 


15± 10 


2 


0.46 


2.2 ± 0.4 


0.33 


M98, IH95, Wa06 




155 


90 


15± 10 


2 


0.58 


2.2 ± 0.4 


0.42 


M98, W03, Wa06 


Hercules 


0.21 




332±io6 


1 




26±^i 




B07, — , S07 


Canes Vn II 


0.071 




336±}i 


1 




38.5ib24 




B07, — , S07 


Leo II 


7.5 


8.7 ±0.9 


35.6±?:? 


1 


0.19 


4.3 ± 2.1 


0.09 


M98* , IH95, K07b 


Leo I 


47.9 


12.6 ± 1.5 


17 ±4 


1 


0.12 


3.4±?:| 


0.07 


M98, IH95, K07a 




76 


13.4 ± 0.7 


5.3 ± 2.6 


1 


0.13 


2.i±l:g 


0.07 


S06, S06, S06 



"Mass-to-light ratios for the Galactic dSphs witii determinations of their internal kinematics. The quoted la error bars do not include the 
associated uncertainties in luminosity or distance. In galaxies in which the reported value of the luminosity varies significantly between 
different authors, we have added a second row. Obviously, in order to have low values of the M/L's, the larger values of the luminosity should 

be taken. 

''We indicate the mass model used. "1" indicates one-single King component models -with or without dark matter. If the velocity dispersion 
profile is available, the assumption that mass-foUows-light can be relaxed and more reliable models can be obtained. "2" refers to models in 
which dark matter and stars follow different profiles (2-component models). 

"^Confidence levels are given in parenthesis. 

* Updated values using D = 233 kpc (Bellaazini et al. 2005). 

Note. — egenda of the acronyms: B06: Bclokurov et al. 2006; B07: Belokurov et al. 2007; IH95: Irwin & Hatzidimitriou 1995; KOI: Kleyna 
et al. 2001; K05: Kleyna et al. 2005; K07a: Koch et al. 2007a; K07b: Koch et al. 2007b; M98: Mateo 1998; M06a: Mufioz et al. 2006a; M06b: 
Munoz et al. 2006b; M07; Martin et al. 2007; OOl: Odcnkirciien et al. 2001; P03: Palma et al. 2003; Q95: Queloz et al. 1995; S06: Sohn 
et al. 2006; S07: Simon & Geha 2007; W03; Walcher et al. 2003; W04; Wilkinson et al. 2004; Wa06: Walker et al. 2006; We06: Westfall et 
al. 2006. 
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Fig. D1. — Comparison of tidal radii in MOND with Ma = 0.7 X lO^"*^ Mq for remote GCs. The error bars enclose la eonfidence limit range. 
The predicted tidal radii (triangles) arc compared with the observed tidal radii (diamonds). The symbols have been shifted horizontally to 
make the plot readable. The observed limiting radii were taken from Harris & van den Bergh (1984) for Pal 14; Harris (1996) for NGC 7006, 
NGC 2419, Eri, and Pal 4; from Hilker (2006) for Pal 3 and AM-1. 
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Fig. D2. — Inferred MONDian mass-to-light ratios versus luminosity for the known Galactic dSph galaxies at 40 < -D < 260 kpc and with 
velocity dispersion estimates (see text). 
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Fig. D3. — The upper panel shows a mass model for Sculptor constrained to agree with observed total stellar content and scale length 
(dotted curve), and a dark matter mass within 1.85 kpc of 5 X IO'^Mq, compatible with the stellar velocity dispersion within this region. The 
core radius of the dark halo is ~ 0.3 kpc. The lower panel gives the contributions to the escape velocity at different radii, from stars (lower 
dotted curve), dark halo (upper dotted curve), and total (solid curve). The dynamics is entirely dominated by the dark matter halo, which 
can only retain LMXBs having a low initial kick. 



